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Abstract 

A symmetry classification of possible interactions in a diatomic 
molecular chain is provided. For nonlinear interactions the group of 
Lie point transformations, leaving the lattice invariant and taking so- 
lutions into solutions, is at most five-dimensional. An example is 
considered in which subgroups of the symmetry group are used to 
reduce the dynamical differential-difference equations to purely differ- 
ence ones. 
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1 Introduction 



The purpose of this article is to analyze possible interactions in a long one- 
dimensional molecule consisting of two types of atoms. The model we con- 
sider is a very general one, described by the equations 

El = Xn- Fniin, t) - GniVn-l, t) = 0, 

(1.1) 

E2 = Vn- Kniin, t) - PniVn, t) = 0, 

where the overdots denote time derivatives and Xn, Vn can be interpreted as 
the displacement of the n-th atom of type X or Y, respectively, from their 
equilibrium positions. We define 

and t is time. The functions F„, G„, Kn and P„ are as yet unspecified 
smooth functions. Indeed, our aim is to classify such systems according to 
the Lie point symmetries that they allow, that is, to classify these functions 
Fn, Gn, Kn and P„. 

The assumptions built into the model are: 

1. The atoms of type X and Y alternate along a fixed uniform one- 
dimensional chain with positions labeled by the integers n 

(see Figure 1). 

2. Only nearest neighbor interactions are considered, i.e. the atom Xn 
interacts only with Yn-i and Y^ and F„ interacts only with X„ and 
Xn+i (see Figure 1). 

3. The system is invariant under a uniform translation of all atoms in the 
molecule and also under a Galilei transformations of the chain. 

4. The systems is strongly coupled, i.e. we assume 

^7^0, ^7^0, ^7^0, ^7^0. 1.3 
din orin-i din dr]n 

5. In the bulk of the article we assume that the interactions are nonlin- 
ear, i.e. at least one of the four functions F„, G„, Kn or P„ depends 
nonlinear ly on the argument ^ or 77, respectively. The linear case will 
be treated separately. 
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6. A discrete symmetry is built into the model. Indeed, the two equations 
are permuted by the transformation 



(1.4) 

p s p r< . ^ J/- . r< p s p 

Models of this type have many applications in classical mechanics, in 
molecular physics, or mathematical biology 0, |^. In applications, the 
form of the functions in eg. (|1 . 1| ) are usually a priori fixed. 

The formalism used in this article is the one called "intrinsic method" in 
earlier articles |§, |^. It has already been applied to monoatomic molecular 
chains p and to a model with two species, or two types of atoms, distributed 
along a double chain [|^. 

In this approach the dependent variables x and y depend on one discrete 
variable n and one continuous variable t. Symmetry transformations, taking 
solutions into solutions, act on the variables x, y and t, not however on the 
lattice variable n. The Lie algebra of the symmetry group is realized by 
vector fields of the form 

X = T{Xn,yn,t)dt + (/>n(x„,?/„,t)<9^,^ + (^^n , 2/n , <9y„ • (1-5) 

The functions r, and ipn are determined from the requirement that the 
second prolongation of the vector field X should annihilate equations ( |1.1|) 
on their solution surface. Explicitly we have ^, |^, |^ 



n+l 

pr(2)X = T{t,Xn,yn)dt+ Yl (pk{t,Xn,yn)d^k 

k=n—l 

n+l 

+ E Mt, Xn, yn)dy, + (j)idi„ + 

k=n—l 



:i.6) 



with 



D^(l)n-{D^T) x„-2(Ar) Xn, 

C = A'^n-(AM?)n-2(Ar)2/„ 

[Dt is the total time derivative). In eq. (|1.6| ) we have spelled out only those 



terms which act on eq.( |l.l|) . 

The use of this formalism is not obligatory. Indeed, the group transfor- 
mations can also act on the lattice ^, 10, 11] and generalized symmetries 
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can be very useful |]T^. In this article we restrict ourselves to the intrinsic 



formalism, described above. 

The present article is organized as follows. In Section 2 we establish the 
general form of the vector fields ( |1.5| ) that realize the symmetry algebra of 
eg . ( |1 . 1|) . We also derive the determining equations for the symmetries and in- 
troduce a "group of allowed transformations" . Allowed transformations take 
equations of the type ( p..l|) into other equations of the same type. They can 
change the functions F„, Gn-, Kn and P„ into other functions of the same ar- 
guments. As in previous articles, we classify equations into symmetry classes 
under the action of allowed transformations 0, |^, 0, |I^ . We also estabhsh 
that equations (|1 . 1|) are invariant under a two-dimensional Abelian group for 
any functions F„, . . . , P„. Section 3 is devoted to Abelian symmetry alge- 
bras. We denote them Aj^k where A means Abelian, j denotes the dimension 
and k = 1,2,3,... enumerates algebras of the same dimension. For each 
interaction we list only the maximal symmetry algebra. Section 4 is devoted 
to nilpotent symmetry algebras, denoted by Nj^^ with the same conventions 
as in Section 3. In Section 5 we find all solvable symmetry algebras with 
non- Abelian nilradicals (SNj^k). In Section 6 those with Abelian nilradicals 
(SAj^k)- AH nonsolvable symmetry algebras are listed in Section 7 (NSj^k)- 
In Sections 3 to 7 we consider only nonlinear interactions. Symmetries of the 
linear case are discussed in Section 8. Conclusions and some applications of 
the symmetries are summed up in the final Section 9. 



2 Determining equations and allowed 
transformations 

The algorithm for finding the symmetry algebra of eq. ( |1.1| ) is 



piX Ea \Et=0 = 0, 



1,2, 6=1,2. 



(2.1) 



The coefficients of all terms of the type ij^ must vanish independently and 
we find that the vector field ( P!75| ) must actually have the form 



X = T{t)dt + 



a + I Xn + Xn[t) 



«+ I yn+ l^n{t) 



yn 1 



(2.2) 
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where a is a constant and A„(t), /i„(t) and r(t) are functions of the indicated 
variables. This form (|2.2| ) is vahd for any interactions F„, G^, and Pn in 
eg. ( |1 . 1| ) ■ Moreover, we have 

T = ro + rl^ + r2^^ (2.3) 

where tq,ti and T2 are constants. 

The constants a, Tj and the functions A„(t) and /i„(t) are subject to two 
further determining equations that involve the interaction functions explic- 
itly. They are 



An + (a - ff) (F„ + Gn) + [A„ - /i„ - (a + f ) F„,5„ 
+ - A„ - (a + 1) G'„,^„_, - T{Fn,t + C^,*) = 0, 

/in + (a - It") (i^n + Pn) + [An - /in " (« + ^n] i^n,^. 
+ [/in - An+1 - (a + i) ^?n] ^n,,;„ " ^(/Cn,* + Pn,t) = 0. 



(2.4) 



(2.5) 



Our task is to perform a complete analysis of eq. (|2.4|) and ( |2.5|) . Concep 



tually, this is very similar to the problem considered in Ref. 7. However, the 
functions figuring in eq . (|1 . 1|) are less general than those of Ref. 7, hence the 
computations are simpler. 

We shall classify the equations of type (|1.1| ) into equivalence classes under 
the action of a group of "allowed transformations". These are transforma- 
tions of the form 

X„ = $n(Xn,ynJ), ?/„ = ^n (^n , J) , t = T{i) (2.6) 

that transform equations (|1 . 1|) into equations of the same form, but do not 
preserve the functions on the right hand side of eq. ( |1.1| ). The requirement 
that no first derivatives should appear and that the transformed functions 
Fn and Kn should depend only on ^n and i, Gn and Pn only on i and r/„_i 
or ?7n, respectively, implies that the transformations actually have the form 

^nit) \ _ -1/2 / Xnii) \ f an(t) A 

= ^^^-7T~^' '^1^4 - C2C3 = 1, q^O, (2.8) 
C3 1 + C4 
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where g, ci, . . . , C4 are constants and a„ and /3„ are arbitrary functions of n 
and t. 

The transformed system is 



with 



(2.9) 



F + (7 

J- n \ 

Kn + Pn 



where 



1,-3/2 
t 



Fn{^n,t) + Gn{r]n-l,t) 



kit) 



, (2.10) 



--1/2 

= yn-Xn = qt {Xn - Vn) + Oi{t) - (3{t) , (2.11) 

r]n = Xn+1 -yn = qt {Xn+1 - Vn) + On+lW " /3n(^), (2.12) 
C4 t — C2 



t 



-C3 t + Ci 

The vector field X of eq. (|2.2|) is transformed into a similar field with 



(2.13) 



r{i) 



T{t{t)) t, a = a, 



U/2 
t 

q 



2' 



an 

fin 



An 



(2.14) 



(2.15) 



The transformed functions and constants must satisfy the same deter- 
mining equations (p.4|) and ( p. 51) . 

As mentioned in the Introduction, translational and Galilei invariance are 
built into the model. That is easy to check. Indeed A„ = ^„ = 1, a = 0, 
r(t) = and = /in = a = 0, rit) = are solutions of eq.(p.4|) and ( p.5| ) 
for Fn, Gn, Kn and P„ arbitrary. No other symmetries exist, unless some 
constraints on the interactions are imposed. 
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We shall use the allowed transformations to simplify the vector fields 
that occur. In particular the coefficient r(t) of a given vector field can be 
transformed into one of the following expressions: 0, 1, tOTf + 1. 

Our strategy will be to first find all Abelian symmetry algebras, then all 
nilpotent (non-Abelian) ones. Once these are known, we can determine all 
solvable ones, having the corresponding Abelian, or nilpotent ones as nil- 
radicals |]15|. Finally, all nonsolvable symmetry algebras will be determined. 



making use of their Levi decomposition [15 



Any symmetry algebra will contain the algebra 

A2,i: X, = d,^ + dy^, X2 = tid,,, + dyJ, (2.16) 
as a subalgebra. Allowed transformations leave the algebra ( |2.16|) invariant. 



Any further element of the symmetry algebra can be transformed into one 
of the following ones 



dx„ + yndy„] 



a = 0, 1, 



Ys 
Y, 



tdt + {a+ -){xnd^^ + yndy^^), 

{t^ + l)dt + {a + t){Xn <9^„ + Vn dyj, 



(2.17) 

(2.18) 

(2.19) 
(2.20) 
(2.21) 



^n{'t){dx„ + dy^), Xn 7^ 0, A„+i 7^ A„, 

Xn(i) dx„ + Xn+lit) dy^, A„ 7^ 0, A„4-i 7^ Xn- 

The interactions that allow these additional terms can easily be determined 
from equations (2^) and (2^). Once this is done, we determine whether the 
considered interactions allows further symmetries. For each interaction, we 
shall only list the maximal symmetry algebra allowed, not lower-dimensional 
subalgebras. 



3 Abelian symmetry algebras 

The lowest dimensional maximal symmetry algebra is A2,i of eq.( p.l6| ), present 
for any interactions in eq. ( |1 . 1| ) . This algebra can be enlarged into a higher 
dimensional Abelian algebra by adding elements of the type ( 2.201 ) or ( p. 21 ). 
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The determining equations for a nonlinear system allow at most four com- 
muting symmetry generators. Moreover, the three-dimensional symmetry 
algebras are never maximal. 

Finally, we obtain two different four- dimensional Abelian symmetry alge- 
bras together with the interactions that allow them. They are 

^4,1 Xi = Ai,„(t)(<9,.„ + 9,j, X2 = A2,„(t)(a.„ + a,j, 

^3 = + dy„ , = t (4„ + <9y J , 

M,n ~ M,n-1 

^^1,71+1 ^l,n 
Al,n 7^ 0, A2,n 7^ 0, Ai,„+i 7^ Ai^ri, 

A2,n _ A2,n ~ A2,n-1 _ A2,n+1 ~ A2,n 
Ai^„ ^l,n ~ Al,„_l Ai^„+i — Xi^n 

X3 = dx„ + dy„ , X4 = t{dx„ + dy,^ ) , 



^l,n+l ~ ^l,n 



K — "^I'^+i f P — P (n t) 

Al,n 7^ 0, A2,n 7^ 0, Ai^„+i 7^ Ai^n, 

A2,n _ A2,n ~ A2,n-1 _ A2,n+1 ~ A2,n 

Ai „ Ai^„ — Ai^„^i Ai^„_|_i — Xi^n 



The algebras ^4,1 and ^44^2 are actually related by the discrete symmetry 
(|1.4| ). Algebra ^4^1 is transformed into ^4^2 by the substitutions 
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Fni^n) > PniVn), Gn{r]n~l) ^ ^n(^n), 

Kn^l{^n-l) ^GniVn-l), Pn-l{Vn-l) > KiCn) , (3.1) 

(^n{t) >(rn+l{t) dy^ , an (t) dy„ >CTn{t) . 

The functions Ai_„(t) and X2,n(t) in algebras A4 i, 2 satisfy the equations 

A2,n _ A2,n ~ ^2,™-! _ A2,n+1 ~ A2,n 

These equations can be solved and we obtain 

Kn = mX2,n + 9{t), A2,„ = - -^d., ^^^^^ 

/(i^) 7^ 0, 7„+i ^ 7„, 

where f{t), g{t) are arbitrary smooth functions of t and 7„ is an arbitrary 
function of n. 

Notice that the quantities Ai,„(t) and A2,n(^) (or /(t), g{t) and 7„) figure 
explicitly in the interaction functions Gn and P„ of ^4^1, or respectively in 
Fn and Kn of A4 2. The two algebras are thus indeed four-dimensional and 
completely specified. 



4 Nilpotent non-Abelian symmetry algebras 

Nilpotent Lie algebras exist for all dimensions dim L > 3. For dim L = 3 
only one type exists, namely the Heisenberg algebra. It has a two-dimensional 
Abelian ideal. Maximality requires that this ideal be the algebra A2^i of 
eq.( p.l6|) . The Heisenberg algebra is obtained by adding the operator T = dt- 
We then calculate the interaction allowing this symmetry algebra, and obtain 

iV3,l Xl=d,„+dy„, X2=t{d,„+dyJ, f = dt, 

Fn = fni^n), Gn = gniVn-l) , 
Kn = ^n('Cn)) Pn = PniVn) ■ 

We mention that this algebra is invariant under the substitution ( |3.1| ). 
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Every nilpotent non-Abelian Lie algebra contains the Heisenberg algebra 
as a subalgebra. We can hence proceed by adding further operators to N^ i. 
Moreover, they can only be added to the Abelian ideal. The determining 
equations (|2.4| ), ( |2.5| ) allow us to add at most two operators. Maximality 
requires that we add precisely two. We thus obtain two mutually isomorphic 
five-dimensional nilpotent Lie algebras with four-dimensional Abelian ideals, 
namely 



XS = + e){d,,^ + dyj, X, = {ant + f + dyj, 

'-'n "^n— 1 

The second algebra N^^2 is obtained from A^s i by the substitution ( p.l|) . 
We mention that the interactions allowing the symmetry algebra A5 i are 
special cases of those allowing the Abelian algebra A^^i. Similarly for A5 2 
and 2. 



5 Solvable nonnilpotent symmetry algebras 
with non-Abelian nilradicals 

A solvable Lie algebra L always has a uniquely defined maximal nilpotent 
ideal, the nilradical NR{L) [|l^]. If a solvable symmetry algebra of the system 
( |1.1| ) has a non-Abelian nilradical, it must be A3 1, A5 1 or A5 2 of Section 4, 
or a four- dimensional subalgebra of A5 1 or A5 2. 

The determining equations ( ^.41 ) and ( |2.5| ) do not allow any extension of 



the four and five-dimensional nilpotent symmetry algebras to solvable ones. 

The Heisenberg algebra N^^i, on the other hand, leads to three different 
four-dimensional solvable symmetry algebras. The Lie algebras are given 
by four basis elements, Ai, A2 and T of A31 and an additional operator Y. 
Below we list these elements Y together with the invariant interactions that 
allow the corresponding symmetry groups. In each case we present a matrix 
A defining the action of Y on the nilradical N^^i. 
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OAT 


Y = 


tdt + (a + ^){Xnd^„ + VnOyJ, 




Fn = 


2a-3 2a-3 




Kn-- 


2a-3 2a-3 

= (^n) kn, Pn = [Vn) ''^+'Pn, 




A = 


diag (a + 2 , 1 , 0- — 2 J5 ^ r "2 ' 2 




Y = 


t9t + (2a;„ + t^)d,^ + {2yn + t^)dy„, 




Fn = 


= fn + l ln(^„), Gn = I ln(r^„„i), 




Kn-- 


= K + ^ In(^n) , Pn = ^ ln(?]n) , 




A = 


/ 2 \ 
12. 


5'iV4,3 


Y = 





A = diag (0,1, -1), cri,„ 7^ a2,n, <7i,n+i 7^ 0-2,n- 

The quantities fn,gn,Pn,kn,cri^n and (J2,n depend on n alone. 

The transformation (|3.1| ) does not lead to any new algebras or interac- 
tions. In the case of the algebra S'A^4,3 we may have a"2,n+i = o"2,n- Then 
(72 can be transformed into a2,n = <J = 0. Similarly, for cr2,n+i 7^ f2,n, but 
c"i,n+i = o"i,n = cr we Can transform into cr = 0. 

6 Solvable nonnilpotent symmetry algebras 
with Abelian nilradicals 

A large number of symmetry algebras of the system ( |1.1] ) is of this type. To 
identify and classify them, we use several known results on the structure of 
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solvable Lie algebras |TH . 



1. The nilradical NR{L) is unique and its dimension satisfies 



dim NR{L) > - dim L. 



(6.1) 



2. Any solvable Lie algebra L can be written as the algebraic sum of 
the nilradical NR{L) and a complementary linear space F, i.e. L = 
F+NR{L). 

3. The derived algebra is contained in its nilradical: [L, L] C NR{L). 

4. For an Abelian nilradical {Xi, . . . , Xn}, the commutation relations can 
be written as 

[-^ii^A;] = iAk)ijXj, [Ak,Ai] = 0, [li, Yfc] = cli^Xi, [Xi,Xk\ = 0, 

(6.2) 

where the elements are the nonnilpotent elements (outside the nil- 
radical). The matrices Ak commute and are linearly nilindependent 
(i.e. no nontrivial linear combination of them is a nilpotent matrix). 
If only one element Y outside the nilradical exists, the nonnilpotent 
matrix A can be taken in Jordan canonical form. 

In our case we can add that the Abelian nilradical must be one of the algebras 
found in Section 3. In principle, the nilradical could be a three-dimensional 
subalgebra of ^4,1, or ^4^2, containing ^42,1 as a subalgebra. However, it 
turns out that all choices of this type lead to symmetry algebras that are not 
maximal for the interactions that they allow. 

The following solvable symmetry algebras occur. 



6.1 dim NR{L) = 2 

The only two-dimensional nilradical that leads to solvable Lie algebras that 
are maximal for the obtained interaction is ^2,1- The solvable Lie algebras 
are always three dimensional. A basis for them consists of Xi and X2 of 
eg. ( |2.16D and an additional element Y, given below. In each case we give the 
matrix A representing the action of Y on the nilradical. 
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SA3,i Xi = + dy„ , X2 = t(a^.„ +dyj, Y = dt + Xnd^^ + yndy„ , 

Fn = Cnfni^n), Gn = Vn-l9n{C,n-l) ^ 
Kn = ^nkriif-^n)-! Pn — VnPniCn) : 

' 1 -1 



Fn = t '^infn{^n)i Gn = t ^?7ri-l(7ri (Cn-l) 5 
-^n i ^n^niy^n) ^ Pn t VriPniCn) i 

UJn = int-^''^'^\ Cn = //nr('^+5), A = diag (a - | , O + i) . 

5^,3 li = a,„ + a,„, X2 = t(a,„ + a,j, y = (t2 + i)ai + (a + t)(x„a,„ + y„a,j, 

F„ = + ly^UnM, Gn = {t^ + l)~'r/n-l/7n(Cn-l), 
Kn = + ly^CnKM, Pn = {t^ + l)-^VnPn{Cn) , 

S^nif' + 1)"^/^ exp[— aarctan(t)], Cn = Vnif' + l)""*^/^ exp[— aarctan(i(:)] 



A 



a -1 
1 a 



These three algebras are nonisomorphic (since the corresponding matrices A 
are not mutually conjugate). Each of these three cases is self conjugate under 
the substitution (R.ll). 



6.2 dim NR{L) = 4 

The nilradical could be three-dimensional, however the obtained solvable 
Lie algebra is never maximal. We only need to deal with four-dimensional 
Abelian ideals of the form A^^i and ^4^2- An extension to a solvable Lie 
algebra is only possible for special cases of the functions Ai,„(t) and X2,n{t) 
figuring in the vector fields and interactions. Below we list all inequivalent 
extensions of A41. There are precisely nine of them. The corresponding 
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extensions of 2 are obtained by the substitution The action of Y on 

{Xi, . . . , X4} is represented by the matrix A. 



VnJ^ 



X3 = cr„e*(a^„ +dyj, Xi = GnQ *(<9^„ + d. 



K„ 



Y = dt + a{xnd^^ + y-nd. 



Vn-1, 



-at 



A 





a-1 








\ 







a + 1 
















a 





\ 








-1 


a J 



SA,,2 Xi = d,„+dy^, X2 = t{d,^+d, 



Vnll 



X3 = cr„ cos(t) + 9y„ ) , Xi = Gn sin(t) + d, 
Y = dt + a{xnd^^ + ?/„5j^„), 



Cnfn (^n 



Pr, 



On ~ On-1 



Vn-1, 



-Or, 



On+1 - On^""' 



-at 



On+1 7^ On 



A 



(a 1 \ 

-1 a 

a 

\ -1 a / 
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Xs = {an + t'){d,^ + X, = {ant + f + 



-^n — ^nkn{^n)i Pn — 





a 








-2\ 




-1 


a 
















a 





I 








-1 


a ) 



y = tc^^ + (a + + Vndy,,), 

Kn^t-^LK{i^n), Pn^a{a-l)t-^ 



Or, 



an+1 -an^""' 



enr^-^+i), an+i^an, a ^0,1 



A = diag (a — q; + |, a + a — i, a+|, a — i). 
y = + (a + + VndyS), 

Fn = t-'Un{cOn), Gn = -\t- \ ^"^ »n-l, 



n = Cnt"^"+5), 7^ an, A 



4 a-„+i - an 

( a -\ 

a 
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\ 




Oa-iy 



a + i 



Xs = a,tV2cos[ln(t)](a,„ + dyj, I4 = a^t^/^ sin[ln(t)](a,„ + dyj, 
Y ^tdt + {a+ + y„dyj, 

-Vn 



(7„+i -(7„''"' 



( a 1 

-la 

a+ I 
\ 



\ 






y = t9t + (a + + ?/„5j/„), 

"-"n '-'n—l 







a-l 





V 








1 \ 
-1 

4-1 

a-i y 
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Xs = Xin{t){dr,^ +dyj, Xin = cTnit^ + 1) exp [« arctan(t)] , 
X4 = X2nit){d^,^ +dy„), \2n = (Jnif + 1) exp [-« arctan(t)] , 
Y = {e + l)dt + (a + + Vndy^), 

F„ = {e + ly^UuM, = (a' + + 77._i, 

= in{t^ + l)"^/^exp[-aarctan(t)], Gn+i 0, 

I a-a 0\ 

a+aOO 

a 1 ■ 
\ -I a I 

Xs = a„(t2 + 1) V2(a,^ + a^J, X4 = an{e + 1)1/2 arctan(t)(9,„ + 9,J, 

F = + 1)5^ + (a + + 

F„ = + l)-2e„/„K), = {e + 1)- ^ Vn-1, 

n '-'n—1 

c^n = + l)"^/^exp[-aarctan(t)], o-„+i 7^ cr„, ay^O, 

/ a \ 

-1 a 

a 1 ■ 

\ -1 a / 

In all cases the interaction terms (?„ and P„ are specified, whereas F„ 
and Kn each involve an arbitrary function of one variable The time 
dependence of the variable Un and the functions F„ and Kn depends on the 
form of the generator Y. 

After the substitution (|3.1| ) we have altogether 18 five-dimensional Lie 
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algebras. No further symmetry generators can be added, at least in the 
nonlinear case studied so far. 



7 Nonsolvable symmetry algebras 

Any finite dimensional Lie algebra L that is not solvable is either semisimple, 
or has a nontrivial and unique Levi decomposition 



L = S>R, (7.1) 

where S is semisimple and R is the radical, i.e. the maximal solvable ideal. 
The only semisimple Lie algebra that can be realized in terms of the vector 
fields ( |2.2| ) is actually simple, namely s/(2,M). Up to allowed transforma- 
tions the realization is unique (and given below by the operators Yi, Y2 and 
Y3). The determining equations ( ^.4| ) and ( |2.5| ) can be used to obtain the 
interaction invariant under the corresponding group SL{2,'R). Equations 



Ll| ) will then be invariant under a five- dimensional group that contains the 



subalgebra ^2,1- We have: 

NS5^i Yi = dt, Y2 = tdt + l{xnd^,^+yndyJ, % = t^dt + t{xnd^,^ + yndyj, 

The Lie algebra NS^^i is isomorphic to the special affine algebra sa//(2,R). 
This is the only maximal nonsolvable symmetry algebra that occurs. 

This completes our analysis of possible symmetries of the system ( |1 . 1|) 
with nonlinear interactions. 

8 Symmetries of linear interactions 

In Sections 3-7 we have excluded the case of linear interactions. Let us turn 
to this case now. We specify equations ( |1.1| ) to be 

Xn = A„(t)e„ + 5„(t)r7„_i + [/„(t), 

(8.1) 
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The system is still strongly coupled, i.e. the functions 5^, C„, D„ are all 
nonzero. The determining equations reduce to 



An - (yUn - An)^n - (A„ - /i„-i)5„ + (a - ^r)Un - rf7„ = 0, (8.2) 

3 

K - (/in - A„)Cn - (A„+i - /i„)D„ + (a - -r)Vn - rVn = 0, (8.3) 

2fA„ + rin = 0, (8.4) 

2fBn + rBn = 0, (8.5) 

2fC„ + rCn = 0, (8.6) 

2fD„ + rZ)„ = 0, (8.7) 

r = ro + rl + r2^^ (8.8) 
since the coefficients of ^n, Vn, Vn-i and 1 vanish separately. 



For An(t), . . . ,Dn{t) generic, we obtain r = and then only equations 
( |8.2| ) and ( |0| ) (with r = 0) survive. These equations can be solved in the 
generic case and we obtain two types of symmetries, both just a consequence 
of linearity. 

1. We take a = and denote {Xh,n, fJ'h,n) the general solution of the ho- 
mogeneous equations, i.e. eq. (|8.1|) with Un = Vn = 0. The vector field 
is 

Xh = Xh,n{t) + AiM(0 dy„- (8.9) 

2. For a 7^ we choose a = —1 and denote some chosen particular solution 
of the inhomogeneous system (|8?T| ) {Xp^n, fJ'p,n) ■ The vector field is 

In particular, if we have Un = Vn = 0, then we take Ap,„ = Atp,n = in 
eq- (|8iq) . 
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The symmetry ( |8.9| ) only means that we can add any solution of the homo- 
geneous equations to a solution of eq.( p.l|) . The symmetry (p.lO|) corresponds 
to the fact that a solution of the homogeneous system can be multiplied by 
a constant. 

Let us now assume that a further symmetry generator exists. It is of the 



form (2.2) with r(t) as in eq.(p.8|). Allowed transformations can be used to 



transform r into one of four cases. Let us consider them separately, 
a) r = 

No symmetries beyond the generic ones are obtained. 



b) r = 1 

Using allowed transformations we simplify the additional vector field into 

f = dt + a{xnd^„+yndyj. (8.11) 

The determining equations restrict the time dependence of the coefficients in 
eq.( ^Tl|) and the system reduces to 

•^n fn ~l~ 9n Vn—1 ~l~ j 

(8.12) 

Vn — ^ra + Pn Vn + ^nC" • 



c) r = t 

The additional vector field and invariant equations are reduced to 

D = tdt + ia+ ^)(x„ + y„ dyj, (8.13) 
" f Q 

k 3 ^^'^'^^ 
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d) r = t2 + 1 

The additional vector field and invariant equations are 

C={t^ + l)dt + {a + t){xnd^„ + VndyJ, 



(8.15) 



Xn = , . N-^ iCn + , ') , x') + — " , , /n exp [a arctan(t)] , 

(r + ir (r + ir {f + lfl^ ^ ^ (8.16) 

i/n = j^r^Y^ in + ^^2^-|^^2 "^n + ^^2 _|!^"]^^3/2 6xp[a arctan(t)] . 

In all cases /„, p^, ti„ and Vn are independent of t. No further 

symmetries exist for any of the interactions ( p.l2[ ), ( ^.14[ ) or ( ^.16[ ). 



9 Conclusions 

Let us sum up the results obtained above. 

For nonlinear interactions the symmetry algebra is at most five-dimensional. 
The following cases occur. 



1. The nonsolvable algebra of Section 7. The dependence of the 
right hand side of eq.( |l.l| ) on and rin is completely specified by 
an inverse cube relation. The dependence on the discrete variable n 
remains arbitrary. The interactions are time independent. 

2. The solvable Lie algebras with Abelian nilradicals SA'^^x^ . . . , 5*^5,9 (and 
SA^^iQ, . . . , 5*745 18 by the substitution (|0| )) of Section 6.2. The inter- 
actions are all "semilinear" . By this we mean that the dependence on 
one variable is specified to be linear, whereas the dependence on 

remains arbitrary (and vice versa for SA^^iq, . . . , 6*^5^18). The time 
dependence of the nonlinear terms in the interaction depends crucially 
on the form of the nonnilpotent element Y. 

Any attempt to enlarge the symmetry algebra by further elements leads 
to linear interactions. 

3. The nilpotent five- dimensional Lie algebras A^s^i and the related alge- 
bra A^5,2 of Section 4. For A^s^i the interaction is again semilinear with 
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Gn and P„ linear in rjn-i and r^^, respectively, and and arbi- 
trary functions of ^„ (and vice versa for N^^2)- The interaction is time 
independent. 

4. Four- dimensional maximal symmetry algebras are either Abelian, or 
solvable with the Heisenberg algebra as a nilradical. For i and 2 
the interaction is semilinear with an arbitrary time dependence in the 
nonlinear terms. For SN^^i, S'A^4,2 and S'A^4,3 the dependence on ^„ and 
rjn is completely specified as being monomial, logarithmic or exponen- 
tial, respectively. There is no time dependence. 

5. A three-dimensional maximal symmetry algebra is either nilpotent, or 
solvable with an Abelian nilradical. For N^^i, the Heisenberg algebra, 
the interaction is time independent, otherwise arbitrary. The model, 
studied by Campa et a/ |]T[], namely 

is of this type. For SA^^i, SA3^2 and 5*^3,3 the interactions involve four 
arbitrary functions of one variable. The interaction is entirely specified 
by the element Y. 

6. As mentioned above, the general interaction in eg. ( P-.l| ) is invariant un- 
der the group of global translations and Galilei transformations, corre- 
sponding to the algebra A2,i of eg. ( |2.16| ). 



The symmetries found in this article can be used to perform symmetry 
reduction on one hand, and to obtain new solutions from known ones, on the 
other. 

Let us look at the example of algebra NS^^i. The system in this 
case reduces to 



^ A _gn_ Yi = — + ^ rg 1) 

'sn 'ln—1 >n In 

The algebra s/(2,M) has three ineguivalent one-dimensional subalgebras, 
namely li, Y2 and ¥3 + l^. Each of them can be used to reduce the sys- 
tem (|9.1| ) to a system of two difference eguations. Let us look at the three 
individual cases separately. 
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Subalgebra Yi 

This algebra leads to stationary solutions. We have 

Xn = Xn,0, Vn = Vnfi (9.2) 

and hence 

f \ 1/3 / h \ 

Jn \ / 



infi = Vn-lfi = Vn,0- (9.3) 

9nJ V 



Subalgebra Y2 

The reduction formulas in this case are 



Xn = XnfiVi, Vn = VnfiVi (9.4) 

and the recursion relations are 



Xn,0 fn 9n Unfi Pn 

^ Sn,0 '/n-1,0 ^ Sn,0 'ln,0 



Subalgebra + 

We put 

= Xn,oVt'^+ 1, yn = ?/n,0 V + 1 (9.6) 

and obtain the recursion relations 

fn 9n kn Pn „\ 

Sn,0 'ln-1,0 Sn,0 '/n,0 

In all three cases we can express ^„ in terms of rjn and obtain a two term 
recursion relation for rjn- These can be solved, but we will not go into the 
details here. 
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■ — — 9-^-= — — ■ — — 9 

Figure 1. Interactions between atoms of type X and Y along a molecular 
chain. 
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